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Abstract: Adapting a method recently proposed by C. Marboe and D. Volin for AA=4 
super-Yang-Mills, we develop an algorithm for a systematic weak coupling expansion of the 
spectrum of anomalous dimensions in the sZ(2)-like sector of planar AA=6 super-Chern-Si- 
mons. The method relies on the Quantum Spectral Curve formulation of the problem and 
the expansion is written in terms of the interpolating function h{X), with coefficients ex¬ 
pressible as combinations of Euler-Zagier sums with alternating signs. We present explicit 
results up to 12 loops (six nontrivial orders) for various twist-1 and twist-2 operators, cor¬ 
responding to triple and double wrapping terms, respectively, which are beyond the reach 
of the Asymptotic Bethe Ansatz as well as Liischer’s corrections. The algorithm works for 
generic states in this sector and in principle can be used to compute arbitrary orders of 
the weak coupling expansion. For the simplest operator with L=1 and spin S'=l, the Fade 
extrapolation of the 12-loop result nicely agrees with the available Thermodynamic Bethe 
Ansatz data in a relatively wide range of values of the coupling. A Mathematica notebook 
with a selection of results is attached. 
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1. Introduction 

The discovery of integrability in the AdS/CFT context [1] has led to important perturbative 
and nonperturbative results for a special selection of quantum gauge models in 2, 3 and 4 
space-time dimensions [2]. The most studied examples are AA=4 super-Yang-Mills (SYM), a 
natural supersymmetric generalization of Quantum Chromodynamics (QCD), and the N=Q 
super-Chern-Simons (SCS) theory in 3d, the so-called ABJM model [3]. Various aspects 
of these systems were scrutinized with the aid of integrable model techniques, such as the 
exact world-sheet and spin-chain S-matrix [4, 5, 6, 7, 8, 9, 10, 11], the Asymptotic and the 
Thermodynamic Bethe Ansatze [1, 12, 8, 13, 14, 15, 16] for anomalous dimensions of local 
gauge-invariant operators, the quark-antiquark potential [17, 18] and the expectation values 
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of Wilson loops in the strong coupling limit [19, 20]. Three-point functions [21, 22] and 
the generalization, beyond the strong coupling limit, of polygonal Wilson loops [23, 24, 25] 
are currently the object of intense research. 

At present, considering the results on perturbative [26, 27], exact [28, 29] and high- 
accuracy numerical methods [30], the spectral problem associated to the study of anomalous 
dimensions in planar AA=4 SYM, can be considered virtually solved. One of the crucial final 
steps was recasting the infinite set of Thermodynamic Bethe Ansatz (TBA) equations into a 
much simpler matrix nonlinear Riemann-Hilbert problem, the so-called Quantum Spectral 
Curve (QSC) [31, 32]. Starting from the TBA equations of [14, 15, 16], this simplification 
was highly nontrivial: it required a deep understanding of the branching properties of the 
TBA solutions [33, 34, 35] and their link with generalized T and Q Baxter’s functions 
[36, 35, 37]. 

Finally, in the paper [26], C. Marboe and D. Volin proposed an iterative procedure for 
the exact determination of perturbative contributions, which was applied to compute up 
to 10 loops for many interesting operators. These results correspond to the sum of a huge 
number of Feynman diagrams and will probably remain - for a very long time - out of 
reach of the more standard approaches to quantum field theory. 

The ABJM theory in 3d is a second interesting example of quantum gauge theory 
which can be considered to be exactly solvable in the planar limit. An exact description 
of the spectrum of anomalous dimensions was obtained by combining information from 
two-loop perturbation theory [38] and on the strong coupling limit, corresponding to the 
classical limit of type IIA superstring theory on AdS^ x CP^ [39, 40, 41]. 

This led to a conjecture for the Asymptotic Bethe Ansatz equations [42], describing op¬ 
erators with large quantum numbers, and to a proposal for the vacuum and the simplest 
excited state TBA equations [43, 44]. The TBA formally encodes the full information on 
the anomalous dimensions in terms of the dressed coupling constant h, the so-called in¬ 
terpolating function which is an essential ingredient of the integrability approach to this 
model. As first noticed in [45], h{X) is a nontrivial function of the t’Hooft coupling A as, 
for example, it scales differently with A at weak coupling {h ~ A) and at strong coupling 
{h ~ y^A/ 2). For the shortest unprotected operator, the so-called 20, the TBA equations 
were solved numerically at intermediate values of h in [46]. 

However, considering the parallel progress made for AA=4 SYM, these remarkable achieve¬ 
ments cannot be considered fully satisfactory and the purpose of [47] was to begin a re¬ 
duction procedure which ultimately led to the Quantum Spectral Curve formulation of 
the spectral problem [48]. Starting from the latter set of equations, exact results on the 
slope function were obtained in [49], together with an important conjecture on the precise 
dependence of h on the t’Hooft coupling A. 

The main purpose of this article is to show that the perturbative scheme described in 
[26] can be adapted to the study of the ABJM theory. We will consider a set of s/(2)-like 
states corresponding to single-trace operators of the form [50] 


tr 




( 1 . 1 ) 


Rigorously speaking, these states do not form a closed subsector but belong to a wider 
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OSp{2\2) subsector. In terms of the global 05p(6|4) symmetry of the theory, they are 
characterised by the Dynkin labels [L + S, S] L, 0, L]; in particular, the operator 20 carries 
the charges L = 1, 5 = 1. The spectrum of s/(2)-like states has been previously studied in 
[13, 51, 52, 53, 54, 55, 56, 46], and explicit results were known up to four loops. 

Although the current work parallels in many ways the N=A SYM case, the difference in the 
analytic properties between the two models makes the study of the AB JM theory nontrivial 
and still very challenging, as both the iterative procedure of [26] and the associated publicly 
available Mathematica program had to be revised step by step. 

The algorithm shows that the answer, at a generic loop order, is expressed as a linear 
combination, with algebraic^ coefficients, of alternating Euler-Zagier sums [57, 58, 59, 60]: 




E n 

l<ni<***<nfc 2=1 


(sgn(aj))’^» 


n. 


\ai\ 


CLj ^ Wj \ {0} , 


( 1 . 2 ) 


where Ofc / 0 to avoid divergence. In contrast, the results of [61, 26] show that anomalous 
dimensions in the sl{2) sector of AA=4 SYM can contain only Euler-Zagier sums with all 
signs strictly positive, or Multiple Zeta Values (MZV). 

Examples of results can be found in Appendix A. We see the appearance of the simplest 
alternating sum, ^_i = — ln(2), at 6 loops for twist-1 operators; irreducible multiple sums 
such as e.g.^ C-i,-3 start to appear from 8 loops, corresponding to double wrapping. 

We point out that the weak coupling expansion of the exact form of h{\) conjectured in 
[49] yields only Q 2 n oc terms (see equation (A.2)); the general transcendental structure 
of the results appears therefore very similar when they are expressed in terms of the t’Hooft 
coupling. Nonetheless, it is certainly worth studying whether this rearrangement of terms 
can reveal some indirect evidence in support of the conjecture of [49]. We leave this for 
future investigations. 

We performed explicit calculations up to 12 loops for the simplest twist-1 and twist-2 
operators, however the procedure is valid for any values of the twist and spin, and can be 
pushed in principle to an arbitrary order. 

Let us collect some comments and observations. We find that the anomalous dimen¬ 
sions do not contain terms of transcendentality one. It would be interesting to see whether 
there are further selection rules for the answers. In the case of N=A SYM, it was observed 
in [61] and [26] that, up to 10 loops, only a subclass of the possible combinations of MZVs 
appear, related to single-valued polylogarithms [63, 64], and it is thus natural to wonder 
whether the same happens for ABJM (either in terms of h or of A). 

Another preliminary observation concerns the shift symmetry: it was noticed in [51] that 
the anomalous dimensions for the operators {L,S) = (l,2n) and {L,S) = (l,2n — 1), 
n G N"'' are the same, at the level of the ABA, up to four loops. Our results suggest that 
this may be the manifestation of an exact symmetry, holding for the maximum and next- 

^In this paper, we restrict to states characterised by a rational Baxter polynomial. In this case, the 
coefficients are actually rational. 

^As a historical remark, this alternating double sum first appeared in the computation of the anomalous 
magnetic moment of the electron at three loops in QED [62]. 
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to-maximum transcendentality part^ of the full anomalous dimension, including wrapping 
corrections. So far we checked this explicitly up to 12 loops for 1 < S' < 6. 

Finally, the analysis of the highest transcendent ality part of the anomalous dimensions 
reveals a significant difference between ABJM and AA=4 SYM. In the latter case, such 
terms are completely determined by the ABA and single-wrapping corrections. This prop¬ 
erty makes it possible to compute this part of the answer as an exact all-loop expression 
[61]. In the ABJM case, we discover that the anomalous dimensions are also affected, at 
all levels of transcendentality, by double and perhaps even higher wrapping corrections^. 
The difference between the two theories in this respect reminds us of the situation with 
the slope function [68], which is a purely ABA quantity in AA=4 SYM but is affected by 
wrapping corrections in the case of ABJM theory [69, 49]. It would be nice to see how the 
computation of [61] could be generalised in ABJM. 

The rest of the paper is organized as follows. In Section 2, we review the QSC con¬ 
struction of [48] and setup some useful notation. In Section 3, we describe the logic of the 
iterative algorithm. In Section 4, we present some checks of the results, against Liischer 
corrections and by comparison with the TBA data of [46]. Section 5 contains the conclu¬ 
sions and some comments on future directions of research. In Appendix A we list the results 
for a selection of operators. Appendices C and D contain technical details completing our 
description of the algorithm. 


2. The Quantum Spectral Curve equations 


In this Section, we will recall the QSC formulation introduced in [48], valid non-perturbatively 
for finite values of the coupling h, which will be the basis of the algorithm presented in 
Section 3. 

The unknown functions appearing in the problem are defined on the complex domain 
of the rapidity u. At finite /i, they in general exhibit an infinite ladder of square-root 
branch cuts, which are described more precisely below. We shall denote by / the second 
sheet evaluation of a function /, obtained by crossing the cut on the real axis, and by /[""I 
the shifted value f{u + obtained by avoiding all cuts. 

The functions {Pi}f=g are particularly simple since, on the principal Riemann sheet, 
they have a single branch cut running along u G (—2/i, 2h). They satisfy the constraint 


(Po)2 = 1 -PiP4 + P2P3. 


( 2 . 1 ) 


The functions instead display infinitely many cuts for u G {—2h,2h) + in, n G Z, 

but have a simple periodicity/anti-periodicity property 


[ 2 ] [ 2 ] 
= CFi^i = , 


( 2 . 2 ) 


®An analogous phenomenon may occur for an asymptotic degeneracy of A/” = 4 SYM considered in 
[65, 66]. The lifting of the degeneracy due to wrapping corrections was calculated in [66], and appears not 
to affect the highest transcendentality part of the result. 

^We are currently unable to disentangle double from higher orders of wrapping. One could try to tackle 
this problem using the generalised Liischer approach described in [67]. 
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where cr = ±1. The necessity to include anti-periodic solutions (with a = —1) was not 
discussed in [48], where the QSC equations were first proposed. Notice however that the 
introduction of this sign is fully consistent with the construction of [48]: in fact, the i^a 
functions were introduced to parametrize quadratically a periodic matrix fj,ab which was 

f 2 l 

proved to satisfy jlab = ^J-ab■ appear very clearly from the weak coupling analysis 

that, in order to capture all states, one needs to allow both signs when taking the square 
root of this relation. The value of a in (2.2) depends on the operator under consideration. 
For instance, for L = 1 one has a = (—1)'^; more generally, a depends on the mode numbers 
of the state and can be defined in terms of the associated Baxter polynomial, as will be 
discussed in Section 3. 

The complete analytic structure is described by a set of Riemann-Hilbert type relations 


Da = -VabX^'^Vc, 


where 


Pab = 


Pab - Pab = 

J^aDb - nDa 


/ 0 -Pi -P2 -Po\ 


1 

0 

0 

0 

Pi 0 -Pq -P3 

oh 

0 0 10 

P2 Pq 0 -P4 

, X = 

0-100 

\Po P3 P4 0 y 


\1 0 0 0 / 


(2.3) 

(2.4) 


(2.5) 


In addition to the previous constraints, we impose some regularity requirements on the 
solutions: the P and the v functions should have no poles on any sheet, and should remain 
bounded at the branch points (locally they behave like \/u ± 2h). 

The global quantum numbers, L G N"*" (twist), S G (spin) and A (conformal 
dimension) for the states described in this paper are encoded in the asymptotic behaviour 
of the solution at large u, through the relations 

Pa ^ AiU'^^, Ao^t°), 

(-A + L- 5)(-A + L + 5-l)(A + L- ,S + l)(A + L + 5) 


^1^14 — — - 


L‘^{2L + l) 


( 2 . 6 ) 

(2.7) 


(-A + L - S' -I- 1)(-A -I- L + 5)(A + L - S'+ 2)(A -I- L + S'+ 1) 

=-(L + 1 ) 2 ( 2 L + 1 )-• 

while the four components of u are required to have the leading power-like asymptotics 


. (2.9) 

The anomalous dimension 7 ~ 0{h?) is defined as in A = L -|- 5 -|- 7 . 

For the following, it is important to underline some further algebraic consequences of the 
equations presented above. First, using (2.3) and (2.1), one finds 

i/a =-Pafe ( 2 . 10 ) 

[ 2 ] 

and, subtracting (2.10), divided by P 12 , to (2.3), shifted by +i and divided by P)^, we 
obtain the following matrix TQ Baxter equation: 

Ti{u)Qb{u) = - (pg(n))”'QPl(n)+ (p'-'](u))"'Q[-2l(n), (2.11) 
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with 


Qa(^x) oc Tl{u) = -a 




(ii) 


pSwJ" ■ 


( 2 . 12 ) 


3. Iteration scheme 

In this Section we describe the iterative procedure for computing the weak coupling ex¬ 
pansion of anomalous dimensions, 


lih) = '^h‘^^‘^^72n+2, 

n=0 


(3.1) 


where, in contrast with the AA=4 SYM setup, each term in (3.1) corresponds to two loop 
orders in ABJM theory (odd loops are vanishing in this model). 

The procedure takes as input the integer quantum numbers L and S. It should be noted 
that there are in general different states associated to a given choice of these charges. Each 
state is expected to be associated unambiguously to a solution of the 2-loop Bethe Ansatz 
[1], and is thus identified by the Baxter polynomial Q{u). As will be explained below, the 
latter is selected at the first iteration of the procedure, from the solutions of the 2 -loop 
Baxter equation (3.20). 


3.1 Analytic assumptions 

The weak coupling expansion is based, as in [26], on a specific ansatz for the form of the P 
functions at finite values of h. It is enough to specify the ansatz for {P}^= 0 ’ since we will 
always consider P 4 to be defined by (2.1). Without loss of generality, we shall set Ai = 1 
and A 2 = h?. This is possible thanks to the symmetries of the QSC equations, analogous to 
the H-symmetry described in [31] (see Appendix B for more details). The choice of scaling 
for A 2 is important, and comes from the observation that, due to the classical value of the 
conformal dimensions. 


A = L + ,S + 0(/i2), 


the combination A 2 A 3 in ( 2 . 8 ) is 0{h?‘) at /i ~ 0 . 
Introducing the Zhukovsky variable 


x{u) 


u -f 'Jv? — 

2h 


(3.2) 


(3.3) 


we can always assume that, on the first Riemann sheet characterised by |a:(tt)| > 1, the P’s 


6 



can be expanded as 


Pi = (xh) ^ 
P2 = (xh)-^ 
Po = {xh)~^ 
P 3 = (xh)-^ 


°° ;,fc\ 

1 + ^ 

V k=l / 

'h “ h>^' 

- + ^C2,k{h) ^ 

V fc=2 


^ 1 

’ 


( ml(u)+ y~]co,fc(/i) 

V A:=l 

/ °° LfeX 

i^2L+l(w) + J^C3,fc(/i) ^ 


(3.4) 

(3.5) 

(3.6) 

(3.7) 


\ fc=i / 

Above, Ml{u) and K 2 l+i{u) are polynomials in u of degree L and 2L + 1 respectively, 
with coefficients depending on h, 


L-l 


2L 


Ml{u) = Ao{h) + '^ mj{h)u^, K2 l+i{u) = A^{h) ^ kj{h) , (3.8) 

i=o j=o 


where Ao(/i) and A^{h) are the same as in (2.6), and we are simply introducing the depen¬ 
dence on h for clarity. It is by determining these coefficients order by order in h that one 
is able to compute the anomalous dimensions using (2.7),(2.8). 

We shall make a crucial assumption on the behaviour of the polynomials, as well as all the 
coefficients of the expansion at weak coupling: 


Ml{u), K 2 l+i{u), Cij{h) 0{h^), for h ~ 0. (3.9) 


This is the same scaling that turned out to be appropriate for the si (2) sector of A7=4 
SYM theory described in [26]. 

To setup an iterative algorithm, we consider the perturbative expansions of the P and 
u functions for /i ~ 0, keeping fixed the value of the rapidity u. These expansions are called 
normal scaling expansions [26] and read more explicitly 


00 


Pz,ns(w) — ^ ^ ^ P2,ns,j'? 

(3.10) 

j=o 


CXD 

(3.11) 


j=0 


where p* is defined by Pj = {hx)~^ Pi. Naturally, the normal scaling expansion of the p’s 
is simply related to (3.4)-(3.7) by expanding at weak coupling both the coefficients and the 
Zhukovsky variable x{u) ~ |(1 -|- 0(/i^)). The procedure also computes the normal scaling 
expansions of p*: 

00 

Pi,ns(«) = Pi,ns,i(w)- (3.12) 

j=0 
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It is important to describe how the analytic properties are modified in the perturbative 
limit. The branch cuts at tt G {—2h,2h) + iZ shrink to zero, and, at every order in the 
expansion, poles normally appear, correspondingly, at positions u G iZ. However, the 
underlying analytic structure of the QSC imposes some very nontrivial constraints. In 
particular, we shall use the fact that at finite coupling the combinations 

Ua{u) +'Pa{u) = Va{u) 

~ f2l 

Uaju) - Ua{u) _ Va{u) - a vj'{u) 

\/v?‘ — y/u^ — 4 / 1 ^ 

are free of cuts and analytic on the whole real axis. This means that the normal scaling re¬ 
expansion of these combinations should not develop a pole around u = 0 at any perturbative 
order. This set of constraints is used in an essential way in the procedure. Besides, with 
the disappearance of the cuts, there is no direct relation between the values of Pj,ns,A:(^) 
and Pi,ns,fc(^) at any given order. Nonetheless, the two expansions are nontrivially related 
through the finite-coupling ansatz (3.4)-(3.7), since, after the replacement x ^ x = 1/x, 
the latter yields an expansion for Pj, which is expected to converge at least in an open 
neighbourhood of |x| = 1 [26]. As we shall see, this will provide precious information for 
the perturbative method. 

3.2 Summary of a single iteration 

Each iteration of the algorithm starts from the knowledge of pi,ns,n and P 2 ,ns,n, comple¬ 
mented with limited information on the normal scaling expansion of pi and p 2 ) and allows 
to compute every other quantity at the same perturbative level. The input needed on pi,ns,n 
and P 2 ,ns,n consists only in the coefficients of their Laurent expansions around u = 0, up to 
the 0(1) or 0(u) term, respectively. For n = 0, all this information can be easily extracted 
from the the ansatz (3.4)-(3.5), and we find 

Pl,ns,0 = 1, P2,ns,0 = 0, (3.15) 

Pl,ns,o(^i) ~ 1 + 0(u), P2,ns,o(^t) ~ « + O(u^). (3.16) 

For the reader’s convenience, we list below the main steps of the algorithm: 

Step 1: Fix the singular parts of po and ps 

At a generic iterative order, the first step is to determine the singular part of po,ns,n and 
P 3 ,ns,n- Since the P’s have a single cut on the principal Riemann sheet, in the normal 
scaling expansion they can develop pole singularities only at u = 0. From the weak 
coupling expansion of the ansatz (3.6), (3.7), we see that, once the pole part is fixed, 
the regular part around u = 0 depends only on the coefficients of the polynomials (3.9). 
This is very important in order to keep the number of unknowns in the algorithm under 
complete control at every iteration. 


(3.13) 

(3.14) 



To fix the pole parts, we consider the sum of equations (2.3) and (2.10) for a = 1, 2, taking 
into account (2.2). Then we obtain 

+ o (po - (hx)^) = p 2 (^ 1^2 + cr 4^'^ - Pi (^I's + a 4^^^ , (3.17) 

(^U 2 + a 4^') (po + (hx)^) = p 3 (^ui + a + pi (^124 + a 4^^) • (3.18) 

Notice that the combinations in brackets are always of the form (3.13), and therefore free 

of poles at tt = 0 at any order. The order-0(/i^) values of the p’s and hx are also regular 

at It = 0. Therefore, the pole part of po,ns,n(^t) and P 3 ,ns,n('w) can be determined without 
requiring the knowledge of t'a,ns,n- 

At leading order, the above computation simply fixes the singular part to be zero. 


Step 2 : Solve the Baxter equation for 124 , with unfixed polynomial part of po ; then fix the 
polynomial part using analyticity 

Then we consider the first of the Baxter equations (2.11), reading explicitly 


,[3] 


,[ 1 ] 



(3.19) 


At leading order in the weak coupling expansion, the rhs can be dropped since P 2 0(^2), 
and (3.19) reduces to the well-known 2-loop Baxter equation ; 


To{u)Q{u) = Q{u + i){u + i/ 2 )^ - Q{u -i){u- i/ 2 )^, (3.20) 

where i^i,ns,o(^i) oc and 

To{u) = a (Mf\u + i/2) - M^^\u - i/2)) , (3.21) 


where M^\u) is the coefficients of hP in the weak coupling expansion of Ml{u). It is 
well-known that the (3.20) admits a discrete set of solutions determining simultaneously 

_ Q 

the transfer matrix eigenvalue Tq{u) and the Baxter polynomial Q{u) = 0^=1 (^ “ ^i)- 
Physical solutions are also required to satisfy the zero-momentum condition (ZMC) [38], 
which can be written as 



kl±|l = ±i, 

Q { 2 ) 


(3.22) 


The sign ambiguity in (3.22) is an important difference as compared to the AA=4 SYM 
case. Notice that the sign on the rhs of (3.22) must be identified with the value of o", 
since requiring the absence of poles in (3.14) at the leading order we obtain precisely 
QW(o)/gf-'](o) = o. In summary, picking a Baxter polynomial corresponds to choosing 
a state in the theory and unambiguously defines the value of a to be used in all subsequent 
iterations. We will discuss shortly how to fix the precise proportionality coefficient a in 
z^i,ns,o(w) = 
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At higher orders, plugging the previously computed values of pi,ns,fc) P2,ns,fc, k < n, 1 ^ 2 ,ns,i, 
I < n — 1, together with the pole part of po,ns,n into (3.19), we are required to solve an 
inhomogeneous Baxter equation of increasing complexity. The technique for its solution, 
coming from [26], is explained in Appendix C. The output of the procedure returns 121 
as a combination of rational functions of u and generalized Hurwitz zeta functions. The 
solution depends on the still undetermined coefficients of the polynomial plus 

an additional set of integration constants Appendix C.3 for details). 

To fix all these parameters, we use the analyticity requirements coming from the regularity 
of (3.13),(3.14) for a = 1, which can be seen as the higher loops analogue of imposing the 
ZMC and polynomiality for Q and Tq in (3.20). This yields a linear system of equations, 
which determines all coefficients, apart from (which is unfixed due to the symmetries 

of the system, see Appendix B) and another coefficient, which is the analogue of a at 
leading order, and corresponds to the freedom to re-define the solution as 1 ^ 1 ,ns,n{u) —^ 
z^i,ns,n + otnQiu). We will show in the next step how to resolve this ambiguity. 

Notice that the size and complexity of this linear system increases with the order of 
iteration. This is one of the delicate parts of the algorithm, since the solvability of the 
linear system requires to use nontrivial relations between the Euler-Zagier sums that are 
the numbers generated by the procedure and make up the matrix of coefficients. Luckily, 
a full set of relations - reducing any multiple sum to a combination of lower weight zetas 
from a minimal basis - have been tabulated up to transcendentality twelve in the datamine 
of [60]®. For the 12-loop computation, we used these results up to weight nine (which also 
critically improves the program’s speed and memory usage); at least one more iteration 
should be relatively simple to perform making full use of these decomposition rules. 


Step 3: Compute p 2 and P 4 ; fix the constants an 
We can now easily determine from the first equation in (2.3); 

= Po zzi - P 2 1/2 + Pi 1^3, (3.23) 


where we know every other term at the relevant perturbative order. The second equation 
from the system (2.4): 

P 2 - P 2 = a . (3-24) 

can then be used to determine p 2 . Considering the small-tt behaviour of P 2 ,ns,n(i^)) and 
matching it with the small-u expansion predicted from the previous iteration (see Step 
6 ), it is possible to fix the coefficient an- In particular, at leading order, the constraint 
(3.16) allows us to fix® 


a = 


h 


-L 


2(gW)2 a„ln|^ 


u=0 


(3.25) 


Finally, we compute p 4 using the quadratic constraint (2.1). Notice that, since P 2 ,ns,o ~ 
this only requires to know ps at the previous perturbative order. 

®In the case of non-alternating sums, these relations are provided up to transcendentality twenty-two. 
®The sign of the square root in (3.25) is unimportant as the QSC is symmetric under Va —>■ —Va- 
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Step 4: Solve the Baxter equation for 1 / 2 , with unfixed coefficients for the polynomial part 
o/pa; impose the analyticity constraints 


We then consider the second of the Baxter equations (2.11): 


; 3 ] 

2 

>[i] 


;-i] 

'2 _ 

>[-i] 


+ a 



[- 1 ]' 

0 _ 

[- 1 ] 


,[1] 


= a 



(3.26) 


Notice that, even at the first iteration, the equation is inhomogeneous. Besides, its ho¬ 
mogeneous part differs from the Baxter equation (3.20) for a sign in front of the transfer 
matrix eigenvalue. 

As explained in Appendix C.2, we can solve this equation and determine 1 ^ 2 , ns, n as a 
function of the yet unfixed coefficients of the polynomial part of ps. By requiring the 
regularity of the combinations (3.13) and (3.14) for 1 ^ 2 , ns, m one again finds a system of 
linear equations whose solution fixes all of the coefficients of the polynomial K 2 L +1 at 
order /i^”, apart from the three coefficients k^^\ which are left free due to 

the symmetries of the system (see Appendix B), and the coefficient of the highest power 


,(2n) 

step. 


The latter contains the anomalous dimension, and will be determined in the next 


Step 5: Compute pi,ras,n and fix the anomalous dimension; compute 124 ^ 

In terms of the previously found solution for 1 / 2 , we can compute pi,ns,n using the first 
equation of the system (2.4): 

Pi - Pi = u . (3.27) 


By matching the expected leading behaviour at u ~ 0, we finally fix the last remaining 
coefficient , and determine the relevant term of the anomalous dimension 72 n-i -2 from 
( 2 . 8 ). 

At this stage we can also compute i^ 4 ,ns,n by inverting the second equation of (2.3): 

a =-P 0122 + ^3 121+Pin4. (3.28) 


Step 6: Compute the seed for the next iteration 

Finally, starting from the knowledge of {Pa,ns,n}J=o’ shall show how to determine the 
quantities needed for the next iteration, namely Pa,ns,n-i-i and the small-u behaviour of 
Pa,ns,n-i-i up to the Same powers as in (3.16) (a =1,2). 

To achieve this, we start by noticing that the ansatz (3.4)-(3.7), which underlines the whole 
method, is naturally organised in terms of the parameter y = h/x. Besides the normal 
scaling expansion, it is very convenient to consider also the double scaling expansion, 
obtained by expanding pi and p 2 at fixed value of y. We can summarise the difference 
between the two expansions by writing 

normal scaling: Pa = Pa 1 . ^ , h], /i ~ 0, u fixed , (3.29) 

\x{u-,h) J 

double scaling: Pa = Pa{y,h), /i ~ 0, y fixed . (3.30) 
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The seed quantities for the n + 1-th iteration will be computed from the truncated double 
scaling series: 


<Sa,n{y) = X] Pa,dsj(y), 0 = 1 , 2 . (3.31) 

j=0 

Let us then show how to obtain (3.31) order by order. The main observation is that, setting 
y —)■ hlx{u) = hx{u) in (3.31) and re-expanding at normal scaling hx{u) ~ u + 0{h‘^), one 
should find precisely the truncated series for pa,ns(^t), up to the same order. The precise 
relation between the two expansions is'^ [26]; 

Pa,ds,niy) = Pa,ns,n('o)|.jj_j^^ — Pa,ns,ni'^)\y^_^y ) O = 1) 2, (3.32) 

where Pansn(u) is the coefficient of /i^” in the normal scaling expansion of Sn-i{hx{u)). 

Relation (3.32) is used in the program to compute Sa^niu) = <Sa,n-i{y) + Pa,ds,n (y)- 
Now, since h/x{u) ~ 0(/i^), it is enough to re-expand Sa,n{h/x{u)) at fixed u to ob¬ 
tain Pa,ns,n-rl(u)- 

Finally, notice that the ansatz (3.4)-(3.7) implies that the small-y behaviour of the double 
scaling expansion is, for every / > 0, 

Pi,ds,z(y) ~ 0{y), P2,ds,z(y) ~ 0(/), (3.33) 

and from (3.32) we see that 

Pl,ns,n+l(u) - p'l,„s,n-rl(u) = 0(u), (3.34) 

P2,ns,n+l(u) - P2,ns,n-rl (^) = 0(u^)> (3-35) 

where Pansn+i(^) be computed from the expansion of Sa,n{hx{u)). This relation 
enables us to fix the small-u behaviour of pi^ns,n-i-i(u), P 2 ,ns,n-i-i(u), up to and including 
the 0(u°), 0 {u) terms, respectively, as needed to restart the algorithm. 

4. Testing the results 

The simplest unprotected operator belonging to the symmetric sector of ABJM is the 
20 with su{ 2 ) and sl{ 2 ) representatives with charges L = 2, S = 1 and L = l,^ = 1 
respectively. By iterating the procedure illustrated in Section 3, we obtained the 12-loop 
result reported in Appendix A. In the same Appendix and in an attached notebook we 
also list the dimensions of a few other operators. 

^Since <S_i = 0, at the first iteration (3.32) is simply 

Pa,ds,o(y) = Pa.ns,o(M)|„^y , 0=1,2. 
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Up to 4 loops, we can directy compare these results with the literature [51, 52, 53, 55, 
56]. In particular, we tested our method for L = 1 and several values of 5 (5 = 1, ...,20) 
against the formula of [51, 56], 

^L=l,s ^ ^ - H. 2,1 + ^-1,-2 - i^-1,2 - i^l,-2 

+ Hi^2 - H 2 -I + i^2,l + ,-1,-1 — ^^-1,-1,1 - ^^1,-1,-1 + .f^l,-l,l^ 

+ Ah^{Hi-H_i)W{l,S) + 0{h^), (4.1) 


where we have used the shorthand notation ilai,...,ai, = ^^ai,...,afc(*S') for the generalized 
harmonic numbers, and the wrapping term yV{L, S) was calculated in [51]; 


W{L,S) = -\ Res 


4^Q 


2 


M=1 


M-l 

X E 

j=0 


q=iM (^2 + M2)iQ ( ^-^(f+1) ) Q Q (^ q+i{M+l) \j 

Q ^ 9-hAf-l)+2ij ^ 




2 


2j — iq — M 2{j + 1) — iq — M 


(4.2) 


The Baxter function to be used in the formula above for L=1 states is 


Q(n) =2Fi(-5,m + l/2;l;2). (4.3) 

The 6-loop contribution, instead, is a new prediction. However, since it is still not affected 
by double-wrapping (that sets in starting from eight loops for twist-1 operators®), we can 
test it against the NLO expansion of Liischer’s formula, which is discussed below. 

In the case of twist-2 operators, we checked our results, up to four loops, against the 
ABA prediction [51, 56]: 

^L= 2 ,S ^ + 8 h^( 2 H _3 + 2 H 3 - H_ 2,-1 - H_ 2,1 - 2//-l,-2 

-2F_i, 2 - 2i?i,_2 - 2Hi^2 - 2H2,-i - i/2,i). (4.4) 

At 6 loops we have the first appearance of single-wrapping effects, calculated in [56] for 
5 = 2,4,..., 30 by using the formula (4.2) with L=2 and 


Q{u) = 3 F 2 {-S/ 2 , iu + 1/2, -iu + 1/2; 1,1; 1), (4.5) 

while the 6-loop ABA result is too long to be written here, but it can be found in [51]. In 
order to check the 8-loop result, we need to expand further the Liischer terms, as for the 
twist-1 6-loop term, involving corrections to the rapidities and the first contribution from 
the transcendental part of the dressing factor: this will be done in Section 4.1. We found 
perfect agreement with all these results. 

®For a state with twist L, single-wrapping terms start to appear perturbatively at order while 

double-wrapping phenomena are expected to contribute from order 
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4.1 More checks with Liischer corrections 

Single-wrapping effects can be calculated at all orders using Liischer’s method [70], first 
generalized in the case of AA=4 SYM in [71, 72], We shall use the following Liischer-like 
formula, introduced for ABJM by [73, 74]: 


^7 


OO 


E 

M=1 


LEM{q) 

27r 




b=l i=l 


(4.6) 


where Saa and Sab are the 517(2j2)-invariant S-matrices found in [75], Ff, = 0,1 for 
bosonic or fermionic indexes, respectively, and is the mirror energy [76, 77] 


Em{q) = 2arcsinh 


[ )■ 


(4.7) 


Using the techniques developed in [72] for the AA=4 SYM case, it is possible to calculate 
the integrand appearing in (4.6) at various orders in the weak coupling expansion. In 
particular, in the symmetric case A = B in (4.6), we can use S-matrix elements squared 
with respect to those used in [72] and the same (fused) scalar factor At the 

leading order at weak coupling, (4.6) reduces to (4.2) for given L and S. More interestingly, 
expanding (4.6) to higher powers of h, we performed further nontrivial checks of the QSC- 
based computation. As discussed in [72], to go beyond the leading order it is essential to 
consider also the finite-size corrections to the physical rapidities or momenta, which yield 
an additional contribution to the anomalous dimension 

5 

hABA = '^iABA{Pk)5Pk- (4.8) 

k=l 


The asymptotic momenta p* appearing in (4.6) and (4.8) are fixed by the Asymptotic Bethe 
Ansatz equations 





s 

n 


X,, 


X, 


x~l-l/x'^x'^ 


(4.9) 


where Q{x^,x^) is the BES dressing phase [8], xf = x{ui Eijl) and p* = —iln^, while 
the corrections (Jp^’s are determined by the system [72] 

s 

L6pk + i'^Spjdp.lnS{pk,Pj) = S^k, (4.10) 

j=i 


where S{pk,Pj) is the scattering phase involved in the corresponding ABA equation, and 
OO « , 4 S 


M=1 ' 


6=1 


1=1 

j¥=k 
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In the case L = 1, 5 = 1, the corrections in (4.8) vanish, since the rapidity is fixed at 
rt = 0 by the trivial Bethe equation = —1, and the contribution to the energy 

determined by the ABA is exactly = \/l + IQh? — 1 at any order. The only 

complication at 6 loops comes from the transcendental part of the dressing factor. The 
latter however can be treated analytically using the techniques explained in [78]. Finally, 
one gets 

= 288 Cs " 768 + 320 C2 + ^ cl + 384 C2 C-I, (4.12) 

while = 256, matching exactly the result reported in Appendix A. 

We also checked the 6-loop term in the L = 1, S' = 2 case, for which we had to calculate 
the corrections (4.8)-(4.11) to the asymptotic momenta 

Pi = -P2 = ^ “ 2/i^ -h I2h^ + 0(/i®). (4.13) 

In particular, we evaluated the integrals appearing in (4.6) and (4.11) using NIntegrate in 
Mathematica, with WorkingPrecision (WP) set to 50, and summed them up to M = 2000, 
obtaining 

76^:'a5n<; = 47.234481782125339488 .... (4.14) 

L=1S=2 

Once the asymptotic energy 7 g = 224 is added, this agrees with the analytic result 

reported in Appendix A up to 18 digits. To improve the accuracy by reducing the errors 
coming from the truncation of the sums, we employed the technique described in Section 
5 of [79]. First, we derived a large-M approximation for the (symmetrised) integrands 
appearing in (4.6) and (4.11). After rescaling q -A sM and keeping the first six orders, one 
finds 

integrandi^ 5 = 2 (s, M) ~ (-1)^ ^ f^A+ 2 L+ 2 n In . (4.15) 

The contribution of the M > 2000 terms can then be estimated by integrating the functions 
an{s) and bn{s) and performing the sum from M = 2001 to oo. Adding this correction, the 
agreement with the analytic result coming from the QSC method reaches 48 digits. Such 
a high precision is sufficient, using EZ-Face [80], to recover the rational coefficients of the 
zeta functions appearing in the QSC result. 

Starting at order /i®, Liischer’s formula (4.6) fails to give the complete answer for 
twist-1 operators, due to the appearance of double-wrapping effects. However, in A7=4 
SYM one has that the highest transcendent ah ty part of the result is affected only by single 
wrapping at any coupling [79, 61]. Based on this expectation, we computed the ABA 
and single-wrapping contribution to the maximally transcendental part of the anomalous 
dimension, at 8 loops and beyond. In the present case, this does not match the result of 
the perturbative solution of the QSC, showing, as anticipated in the Introduction, that 
these terms also include double wrapping corrections. To gain more insights, it would be 
interesting to extend the Liischer-like techniques to the double-wrapping order, using the 
method developed in [67]. 
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For twist-2 operators, instead, double wrapping effects start only at 12 loops, and we 
can still test the 8- and 10-loop results. In particular, for 5=2 we obtained 

TT /i2 23/1^ 8/i‘^(20-b9C3) /i8(9857-b 5760 C3 +8640Cs) 1n^ 

” = + ^-+-Sw!-+ °''' >• 

(4.16) 

from the perturbative expansion of the Bethe equations (4.9). Then the ABA contribution 
to the energy is 

2 

'JABA = \/l + sin(pfc)2 - 1 (4.17) 

k=l 

= Ah^ - 8/i^ + 32/i® - 4(39 + 8 C^)h^ + 4(209 + 72 C 3 + 80 + 0(/i^2). 

Plugging the solutions (4.16) into (4.6) and expanding at leading order, we get the expected 
result 

r —9 Q—o 16 64 

76 = -y + y C 2 - 42 C 4 . (4.18) 

At 8 -loop, evaluating numerically with WP=50 and summing up to M = 2000 the integrals 
coming from the expansion of (4.6) and (4.11), we obtain 

IsZlfppLg = 302.916652975526460098546231..., (4.19) 

that, added to the corresponding terms of (4.17), agrees with the result reported in Ap¬ 
pendix A up to 25 digits. A 49 digits agreement is obtained by adopting the asymptotic 
resummation method described above. 

Since the complexity and the size of the integrands involved in these calculations 
increase drastically with the value of 5 and the loop order, we managed to perform just 
some lower-accuracy checks for the L = 2, 5 = 2 operator at 10 loops and the L = 1, 
5 = 3,4 operators at 6 loops. Keeping the first 200 terms in the sums, we found 

TfoSSna = -4711.383655826655401..., = -35.71109708632629843 ..., 

(4.20) 

which both agree with the results reported in Appendix A up to 17 digits, once we 
add 7 iq~Jba^ = 4(209 -|- 72^3 -|- 8 OC 5 ) and = 82720/243 respectively. Finally, 

including the first 50 integrals for the case L = 1, 5 = 4, we obtained = 

343.369048487641320..., matching the expected result with 16 digits precision. 

4.2 Comparison with TBA 

The exact anomalous dimension of the 20 operator at finite values of h was studied by 
F. Levkovich-Maslyuk in [46], up to values of /i ~ 1, by solving numerically the TBA 
equations proposed in [44]. In Figure 1 is shown the interpolation of the TBA data of [46], 
together with various truncations of the perturbative result, 

7^=1’'^=!(h) = 8 / 1 ^ - 26.3189-h 195.487 h® - 1966.44/i® -h 22309.5- 271422. + ... 
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y(h) - yABA(h) 



Figure 1: Comparison of perturbative results with the TEA data of [46]. Notice the excellent 
agreement of the non-perturbative TEA result with the [6/6] Fade approximant curve obtained 
from the 12-loop dimension. 


The plots show the difference between the anomalous dimension and the exact ABA result 
= vT+TfT/? — 1. An order [6/6] Bade extrapolation of the 12-loop result yields 


Fade 


7 


L=1,S=1 




38.8151/i^ -b 5.6811 


164.024 + 81.7529 + 17.4840 h? + 1. 


, (4.21) 


which is in remarkably good agreement with the TBA data up to /i ~ 0.4, as can be seen in 
the figure. This goes beyond the radius of convergence of the perturbative series, which - 
by analyticity arguments supported by a fair amount of numerical evidence in A/=4 SYM 
[8, 81, 26] - is expected to be j/ij < j. The approximation is still rather good decreasing the 
number of loops, e.g., with the optimal choice of the order of the Fade approximant, the 
relative error in the estimate of \'^{h) — 7aba(^)|, compared to the TBA value at /i ~ 0.3, 
is 2.4 X 10“^, 1.2 X 10“^, 2.1 x 10“^ for the Fade extrapolation of the 8-, 10-, and 12-loop 
results, respectively. 


5. Conclusions 

The discovery of integrability in the context of the AdS/CFT correspondence has opened 
the exciting opportunity to solve certain interacting gauge theories in 3d and 4d. Within the 
powerful integrability setup, exact expressions for particular physical observables are often 
fairly easy to obtain, revealing the beautiful simplicity of the final outcomes. One of the 
recent achievements in this research field is related to the study of anomalous dimensions 
and their weak coupling expansion through the so-called Quantum Spectral Curve. 

In this article, relying on the success obtained for A/'=4 SYM, we have developed an 
algorithm to solve perturbatively the QSC for the AB JM model, finding a perfect agreement 
both with the TBA numerical results and with earlier perturbative computations and 
predicting several new terms. We are planning to optimise the Mathematica code and to 
render it publicly available. 
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Along the way, we have recorded some important formal differences between the present 
case and AA=4 SYM. These include the appearance, for the ABJM model, of Euler-Zagier 
sums with both positive and negative signs, MZVs of even arguments C 2 n and the experi¬ 
mental observation that the highest transcendentality part of the anomalous dimensions is 
not completely determined by the ABA and single wrapping corrections but also contains 
double (and possibly higher) wrapping contributions. 

There are a few interesting open problems related to the current work. First, as already 
mentioned in the Introduction, it would be interesting to understand better the structure 
of the perturbative outcomes, and investigate whether some patterns emerge. One could 
also try to use our results to deduce expressions for the anomalous dimensions for all values 
of the spin, as was done in [27] in the case of AA = 4 SYM. This would be particularly 
interesting for the purpose of exploring various conjectures on the large-spin behaviour 
[51, 82]. The generalization of this perturbative approach to other sectors of the theory 
is another obviously important problem. Although we do not have yet concrete results 
to discuss, the generalization to the full OSp{2\2) sector should be fairly straightforward, 
while setting up the iterative procedure for generic operators appears to be more involved. 
Finally, it would be very nice to complement the current results by adapting to ABJM the 
numerical technique developed in [30] for the non-perturbative solution of the QSC, or to 
try to transfer some of these powerful methods to the study of the QSC for the Hubbard 
model [83]. 
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A. Sample results 

In this Appendix we present some explicit results, which can also be found (including three 
more operators) in the Mathematica notebook Results.nb attached to the present paper. 

Using reduction formulas such as the ones in [84], the result can be written in terms of a 
small number of non-reducible sums. We will use the basis of [59, 84], which is conjectured 
to be minimal. We find empirically that the /-loop result involves only the basis elements 
Cai,...,afe with weight w = Yli=i l®*l < / — 2 and depth k < I — I — w. In particular, the 
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12-loop result can be written in terms of the following sums: 


C-l) C 2 , {C2n+l}^=l , 

Cl-3, Cl-5, C 3 - 5 , Cl-7, 

Cl,1,-3, Cl,1,-5, Cl,3,-3, 

Cl,1,1,-3, 


We report the results in terms of the interpolating function h{X). Assuming the validity 
of the conjecture of [49] for this quantity, namely, 


, sinh(27r/i(A)) ^ /I 1 1 3 , 0 . , , 

A =--3O I 2 . 5, 5; 1 . 5; - snih (2»/,(A)) I . 


(A.l) 


the anomalous dimensions can be rewritten in terms of the true coupling constant A using 
the expansion: 


= A (1 - 2 (A^,) + 15 (A^ C3)^ - (A= 6)’ + (A^ 6)‘ - (A^ Ci) 


50501 


88949769 


2 A a5 


+ ..., 


(A.2) 


the first two orders of which have been verified by direct calculations [52], 


Twist-1 

For the 20 operator with L=l, 5=1, Q{u) = u, we found: 

7L=i,s=i = h^S-h^l 6 C 2 + h^ f + 288Cs + 384C 2 C-i + 320C 2 - 512 J 

+ (rs=i,2 - - 4736C3 - 1536 C 2 C-i - 512 C 2 + 8192) 

+ + 13632 Cs + 60928 C 2 Ca + Cl C-i + 12288 Ci,-3 + 21504C 3 C-i 

-t 35840 Cl + 6144 C 2 C-i - 61440 C 2 C-i - 25600 C 3 - 102400 C 2 - 32768) 

1,12 Ct3(i2) 559684032CI 11225408(7 1509696(1(3 25762912(2(5 2770944 ^ 3 ^ 

+ * + u + u ^C2(-l 

39740672 

- 245760 (2 ( 1,-3 + 122880 (i ,_5 - 682496 (| - 155136 (5 (-1 - 1523712 (2 (3 (-1 - 

lUo 

1701 

-^-(I (-1 - 98304 ((i,_3 (-1 - (i,i.-3) + 145920(5 - 1164288(2 (3 - 86016(3 (li 

- 16384(2 (ii - (I (_i + 688128 (i,_3 + 1204224(3 (_i + 

-t 344064(2 (li -t 2441216(3 + 1769472(2 (-1 -t 2490368(2 - 1572864) -h 0(h^^), 
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and for L=l, 5=2, Q{u) = — 4: 


7 


L=1,S=2 _ ,2 


32 


= h^8-h^ ( —+ I 6 C 2 ) 

O 


5 3 3 


+ h^lT 


8 f -p(8) 


21184 Cl 


S=l,2 

lIO fp(10) 

“ S=l,2 

380224 Cl 
4^ 


15 

351264 Cl 
35 


- 1792C3 + 


4352 C 2 29696 \ 


27 / 

172032 


5952 C 5 + 38016 C 2 Cs + Cl C-i + 4096 Ci,-3 + 7168 C 3 C-i 


+ 2048 C 2 Cli - 22528 C 2 C-i - 


128000 Cs 270848 C 2 55808 \ 


{ p(12) 

S=l,2 


9 9 27 J 

917102848 Cl 8888248 Ct 9522816 CICs 20301152 C 2 C 5 6177792 ^3 
9975 17 85 17 35 


CIC-I 


401408 C 2 C 1.-3 286720 C 1.-5 705280 C| 1735168 1495040 

- O - o - Cs C -1 -^- C 2 Cs C -1 


- 77^ Cl <i. + 65536 L,_3 - d C-.) - 


- 57344C3 Ci, - 7^ C 2 Cl, + Cl C-1 + 212992 Ci,-3 + 372736C3 C-i + 


1760896 Cs 1138432 C 2 C 3 
9 27 

8532992 C| 


106496 C 2 Cli + 278528 C 2 C-i + 


15 

15569920 Cs 11159552 C 2 54956032\ 


45 


27 


27 


243 / 


0 {h 


14 N 


where for brevity we have grouped together a number of terms that are common to the 
5=1 and 5=2 cases, at orders /i®, and /i^^; 


,( 8 ) 24672 Cl 

35 

,( 10 ) _ 911816 Cl 


- 3520 Cs - 1920Cl C -1 + 864C 2 Ca - 3072 C 1.-3 - 5376C3 C-i - 1536 C 2 Cli, 
78336 


S=l,2 


175 


34860 Ct 


Cl C -1 - 7120 C 2 Cs - 128 Cl Ca + 10240 Ci.-s 


- 1024 C 2 Ci.-a - 992 C| + 19840 Cs C-i + 32768 C 2 Ca C-i + 19968 C| Cli + 21504 Cs Cli 

- 24576 (Cm,-a - Ci,-a C-i) + 4096 C 2 Cli, 


11002416 Cl 


o 125088 Cl Ca 80792 CiCs , 

2 - 335664 Co--+ 67340 C 2 Ct -^ Cl C-i 


918336 


35 

10240 C 2 C 1 .-S 


5 7 

3358720 Cs,-s 2099680C3 Cs 
m 19 


p(12)_ 

“ 275 

12707840 Ci,-T 29696 C| Ci.-a 
57 5 

+ 124672C 2 Cl + 142240Ct C-i - 496640C 2 Cs C-i - Cl Cs C-i - Cl Cli 

1589248Cl,s.-s /■ 3 - , 5357568Ci,i.-s 

-- 235520 Cl ,-3 Ca ~ 565248 C 2 Ci,i ,-3 + 245760 Ci ,-5 C-i 4- ^ - 

+ 238080 Cs Cli - 688128 C 2 Ca Cli - 636928 C| C-i - 135168 C| Cli - 98304 Cli Ci.-a 
+ 196608 (C -1 Ci,i.-a - Ci,i,i.-a) - 8192 Cli C 2 - 57344 Cli Ca- 


We shall give the following results up to double wrapping. For L=1 and 5=3, Q = u^ — \u: 


L=1,S=3 ^ _ ^4 


7 

+ h^ 


672 Cl 


78) 

5=3,4 


3 

+ 512C3 


224 64 C 2 

2048 , , 

— 5 — C 2 C -1 + 


48448 C 2 590944 \ 


135 


1215 ) 


2647232 Cl 2658304 Cs 45056 , . 10672576 C 2 1942407808\ 

' 52C -1 + 


675 


405 


135 


6075 


273375 I 


0 {h 


ION 
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and for L=1 and 5'=4, Q = 


L=1.S=4 ^ ^ . 4 

3 




V 


512 C3 


p872 64(2 

2048 , , 

C2 C-1 + 


151744(2 1033504 \ 


3 ’ ’ 945 8505 

7831552(3 2428928(2 (-1 781130944(2 


2835 945 

where we have used the shorthand 


297675 


+ 


( _ 

S=3,4 

47947882624\ 
93767625 / 


13464896 

4725 

0(0, 


^,8, 32896CI 56320<b , , 10240 ,2, 65536 , 

0^3,4 =-35-g— + 1536C2C3 - ^ C2 C-. - ^ C.,-3 

114688 , , 32768 ^ ,2 

-^C3C-i-^C2(!i. 

Up to the loop order we have reached, these dimensions evaluate numerically to 


^L- 1 , 5-1 2 ;. 8/i^ - 26.3189/1'^ + 195.487/4® - 1966.44/4* + 22309.5/4^° - 271422./4^^, 

^L=i,5=2 _ 30 9850 /^4 271.234 /4® - 2562.52 /4* + 28298.5 - 340713. /4^^, 

32 

^L=i,S=3 ^ _ 40,0697 /4^ + 304.700 /4® - 3019.43 /4* + 33855.8 - 408283. 

32 

^L=i,5=4 235 _ /j2 _ 47,3840 + 343.369 /4® - 3266.52 /4* + 36262.1 - 436578. 

More precise numerical values can be found in the attached notebook Results.nb. 


Twist-2 


For the simplest twist-2 operator with L=2 and S=2, Q = v? — a = -|-1: 


7^=2.S=2^4^2_8^4_;^i 

^/^s (1240(1 
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which evaluates to 


^L=2,s=2 4 /i 2 _ g^4 16.301018/1® + 108.45083/1® - 3197.3744/i^° + 56169.236/i^^. 


B. Symmetry 

In this Appendix we shall discuss a simple symmetry of the QSC equations, which allows 
us to choose freely the constants Ai and A 2 and four of the coefficients of the polynomials 
(3.8). The simplest symmetry of the QSC equations (see [32] and Section 2.3 of [28] for 
the very similar AA=4 case) is the transformation 

I/b, Pab —t Ra Pjj R^, (BT) 

where R is any 4x4 constant matrix satisfying 

R^xR = X- (B.2) 


To preserve the structure of the QSC, we should also require that the transformation does 
not change the ordering of the magnitudes of the P functions at large tt, jP 2 | < |Pi| < 
|Po| < IP 4 I < [Psl- The most general form of R compatible with these constraints has six 
degrees of freedom and can be written as: 


R = 


ai 

0 

0 

0 

a2 

as 

a4 

0 

_25 

0 


0 

«3 

«5 

«3 

02 + 040.5 

j_ 

Oiias 

Ql 

aias 

01 


and the transformation acts on the P functions as 

Pi —^ ai (asPi + a4P2), P 2 —^ — P 2 , 

as 

P . p I n P I “2 + «4a5 p 

■VQ —^ -VQ + as Pi H-P 2 , 

as 

as , «4 ^ , o “2 ^ ^ ^ , «2 - Oiia^aQ 

Ps — Ps H-P 4 + 2 — Po - ae Pi H-P 2 , 

ai ai ai aias 


1 


ai as 


(P 4 — 2 as Pq — as Pi — (0405 + 2 a 2 as + aiag) P 2 ) • 


(B.3) 


(B.4) 

(B.5) 

(B. 6 ) 

(B.7) 


From (B.4), it is simple to see that the constants Ai and A 2 characterising the leading 
asymptotics Pi ~ Aiu~^, P 2 ~ A 2 U~^~^ can be fixed to arbitrary values by an appropriate 
choice of ai and as- Similarly, relations (B.5),(B. 6 ) show that the coefficients mo{h), ko{h), 
kL{h) and A: 2 l(/i) defined in (3.8) - which correspond to certain coefficients of the large-rt 
expansion of Pq and Ps - can be chosen freely by tuning as, oq, a 2 and a 4 , respectively. 
Accordingly, these numbers are not fixed by the algorithm at any order in h. 
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C. Solving inhomogeneous Baxter equations 


In this Appendix, we shall present the basic method to solve the inhomogeneous Baxter 
equations (3.19) and (3.26) encountered in the procedure. 


C.l Solving the inhomogeneous Baxter equation for I'l 

At a generic perturbative order 0{h‘^^), equation (3.19) reduces to 

(u + i/2)^ -{u- i/2)^ -Toqi = (C.l) 

where q^ = 17i a source term of increasing complexity, and Tq is the zero-order 

transfer matrix. At leading order, we know that the source term is zero, and the regular 
solution is the Baxter polynomial Q. To solve the generic case, we follow the method of 
[26]. Considering the ansatz qi = Q it is simple to see that, in order for (C.l) to be 
fulfilled, / must satisfy 

V_ (^r^QWg[-ilv+(/i)) = Ui qW, (C.2) 

where we have denoted V+ g = g — g^^^ and V_ g = g + We introduce the inverse 
operators and T-, such that = g. We shall give a precise operative definition 

of the operators in Section D.2, by explaining how they act on the functions generated 
by the algorithm. Using these operators, an inhomogeneous solution of (C.l) can be found 
as 

/l,inhomo = «'+ (^Lg[i]g[-1]^- • (C-3) 

[21 

Moreover, setting Ui = 0 and 'I'_(0) = <ki^anti, where <l>i^anti = “*^1 anti ^ generic anti¬ 
symmetric function, we find a second independent solution of the homogeneous equation, 


/ijhomo — ^l,anti 'k— 


Introducing for convenience the notation 




^[“^1 = g[“^i dz ( _^^ 


and putting all pieces together, the general solution of (C.l) can be written as 

1 


= ^l,perg'"^' + ^2,anti ( 


n^gwg[-ii 


Ui 


(C.4) 


(C.5) 


(C.6) 


f2l 

where 4>i^per = denotes a generic Aperiodic function. Following [26], it is convenient 

to rewrite this expression in order to cancel its apparent poles at the Bethe roots where 
g = 0. To achieve this, we introduce two polynomials A and B, of degree 5 — 1, defined 
by the condition 


AgW + Bg[-^i = 1. 


(C.7) 
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The Baxter equation then implies that 


= QR, 


(C.8) 


with R a polynomial of degree L — 2. Introducing the constants r±^k and the polynomial 
C through 


R 


E 


A 




rk,+ rk- \ 

(udl)^ (u[“b)fcy ’ 

(C.9) 


(C.lO) 


k=l 


we see that (C.5) can be written as 


- 2 ^^ + + rk-)i]_k(u)j , 


(C.ll) 


where ri_k{u), k G N"*" are certain polygamma functions (see equation (D.4) below), solu¬ 
tions of 


, [ 2 ] 1 
V-k + rj-k = -k- 


(C.12) 


Likewise, the inhomogeneous solution of (C.3) can be written in the following form, which 
only involves poles at positions u G iZ: 

/l,inhomo = (c/igW) C + Q[-d4/+ (c/iQW) ^ • 


(C.13) 


C.2 Solving the inhomogeneous Baxter equation for U 2 

The equation arising from (3.26) at the (n -|- l)-th iteration of the algorithm has the form 

{u + i/2)^ qf -{u- i/2)^ 4“'' + To 92 = (C.14) 

with 92 = 44s n- cine can proceed in a similar way as for (C.l), by paying attention to the 
different sign in front of the Tq 92 term. This implies that a simple, homogeneous solution 

f 2 l 

of the equation for t /2 = 0 is simply given by 92 = g<h 2 ,anti, where ‘h 2 ,anti = 

any anti-periodic function. Likewise, we see that an independent family of solutions to the 

homogeneous equation is described by 

92 = 44ir^> (C-15) 

\ 2 \ 

where 4 > 2 ,per = *^2 per ^-periodic and Z is defined in (C.ll). 

Making the ansatz 92 = Q /4 > di^d 

v+ (u^gWg'-^'v_(/2)) = -i/sgw, (c.ie) 
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so that a solution to the inhomogeneous equation can be found as. q 2 = Q f^inhomo^ with 


/2,mhomo - , (C.17) 

and the most general solution can be written as 

4"'’ = + <^2,perZ + '&+ ( ^Lq[1]q[-1] 

More explicitly, in terms of the same polynomial C and constants rk,± defined in (C.9), 
(C.17) can be expressed as 

/2,inhomo (C.19) 

= (-t/sQW) C + 1^^/+ (-t/sQW) ^ - C[21c/2^ . 

C.3 Periodic coefficient functions 

To construct periodic/anti-periodic functions without introducing unphysical poles or an 
unphysical exponential growth at infinity, we consider the following combinations 

Vk{u)=i^k{u)+sgn{k)4k'^\u)=sgn{k)Vk{u + i), 0 / A: G Z, (C.20) 


where r]k is defined in (D.3). The coefficient functions appearing in (C.6), (C.18) are then 
constructed, at every iteration of the program, as 




A 


(C.21) 


for a = 1,2, where parameters. The number of terms included in the 

sum increases with the perturbative order: at the n-th iteration, we may take A = 2(n — 1). 


D. Functions generated by the algorithm 

In this Appendix, we give a precise definition of the operations This will show that 
the algorithm always produces answers in a specific algebra of functions comprising: 

• rational functions of u, with poles only in the set u G iZ ; 

• the functions tja, with A a multi-index A = oi, 02 ,..., a^, a, G Z \ {0}, which are a 
generalisation of Hurwitz zeta functions [61] , 

• periodic/anti-periodic functions constructed as explained in Section C.3. 

The presence of alternating signs is an important difference as compared to the A7=4 
SYM case, leading to the appearance of alternating Euler-Zagier sums in the results. 
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D.l The 77 functions 

When the sums are convergent, the operators and can be implemented as 

00 00 

^^{g) = Y, ^-( 5 ) = (D.l) 

n=0 n=0 

This leads to the natural definition of the r/a functions satisfying V±? 7 ±|(i|(ii) = 


00 

Vaiu) = E 

n=0 


(sgn(a))^ 
(u + in)l“l ’ 


which is simply expressible in terms of polygammas: 


(D.2) 


^a(u) 


V-a(u) 


(a-1)1 




a G N+, 


2 “ (a - 1 )! 







a G N+. 


(D.3) 

(D.4) 


It should be noted that the sum (D.2) is divergent for a = 1, but we may use (D.3) as 
regularised definition of r/i. The ambiguity in this choice does not affect any physical 
results. In fact, notice that the g functions enter the algorithm through the solution of 
inhomogeneous Baxter equations and this ambiguity only goes into a redefinition of the 
integration constants described in Section C. Another important point to underline is that 
the particular regularisation of gi introduces igi (i) = = 7Euier-Mascheroni in the solution 

of the QSC. While this number enters the expansions of various P and u functions (which 
are not directly physical due to the symmetry described in Appendix B), it always cancels 
out of the anomalous dimensions. 

We then define ga,B, with B any multi-index, as follows 


Va,B{u) = E(sgn(a))^ ^^ , , . 

n=0 ^ ' 


Explicitly, these functions can be written as 


'nai,a2,...,akiu) — E! 0^ 

nk>nk-i>--->ni>0 i=l 


(sgn(ai))”» 


ni—ni—-\ —1 


{u + 


(D.5) 


(D. 6 ) 


where by convention we set no = —1. The sum is convergent for 7 ^ 1. One can define 
the marginally divergent cases where = 1 through the stujfle algebra (see Appendix A 
of [61]), which allows one to express them in terms of convergent g functions and gi. 

Relation (D. 6 ), shows that the Laurent expansion of the g functions around u = i 
naturally produces multiple Euler-Zagier sums, which therefore enter the algorithm. The 
precise relation is 

gA{i) = i-^^^NACA, (D-7) 
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where, given the multi-index A = {ai,..., Uk}, we have defined 


k k 

\^\ = ^Wk\, = ]Jsgn(aj), (D.8) 

t=i i=i 

and the multi-index A is 

A = {ai sgn(a2) , . . . , Qjyyi sgn{am+i),.. ■ ,ak} ■ (D.9) 

Alternatively, it is useful to note the following relation with the harmonic polylogarithms 
Ha{x) of [85] evaluated at unity (see also Section 10 of [86]): 

= (D.io) 

where A is obtained by reversing the order of the indices in A. Relation (D.IO) is the most 
useful for comparisons with the notation of the datamine [60]. 


D.2 Defining the operators dl'j- 

Let us finally come to the full definition of the linear operators Applied to a polynomial, 
we require that they yield a polynomial answer of the form 




E 

vi=o 


a,' = 




tti = 


i=o 




m+1 

E 

i=i 


Co w. 


(D.ll) 


Each rational function of u is then broken into a polynomial part and a sum of inverse 
powers, leading to the appearance of r] functions 




1 

mI“I 


ri±\a\- 


(D.12) 


As in [26], we then define recursively the action of on the products of rational functions 
and ? 7 ’s. One starts form noticing the following relation (which is a consequence of the 
nested definition (D.5)) 


Va,A 


[ 2 ] 

(sgn(a))r?lX + ^, 


for any multi-index A, which leads immediately to 




nr-' \ 

{u + i n)l“l j 


[2n] 

^±\alA- 


(D.13) 


(D.14) 


Besides, when more general combinations like /{u+i /c)l“l are encountered, one can use 
(D.13) until either a term of the form (D.14) is met, or the r] runs out of indices ( 7/0 = 1). 
To deal with the polynomial-times-?? part, we use the relation 




uH??H±(sgn(6))^±ignw(^'“')''’ 


„[2+2n] 

VA 


(u -I- in)\^\ ’ 


(D.15) 
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which can be checked by a little algebra, and finally obtain 



=F (Sgn(6))^'± 


„[2+2n] ' 

{u + in)\^\ 


(D.16) 


The i-periodic/antiperiodic functions Vk are easily dealt with using 

'^±{'Pk ■■■) ='Pk T±(ggn(k)) (...)• (D.17) 

In conclusion, the operators are fully defined on the algebra of trilinear combinations 
of rational, rj and Vk functions. We notice that, since some of the equations adopted in the 
algorithm are quadratic, products of two rj functions may also be generated along the way. 
However, by using stuffle algebra relations such as the ones in [61], these can be converted 
into expressions that are linear in all the rj functions. Up to the order we reached, it actually 
turned out that this step was unnecessary, since all the source terms in the inhomogeneous 
Baxter equations were already linear in the ry’s. 
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